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Abstract 

We study the dynamics of the linear and non-linear serial dependencies in fi- 
nancial time series in a rolling window framework. In particular, we focus on 
the detection of episodes of statistically significant two- and three-point cor- 
relations in the returns of several leading currency exchange rates that could 
offer some potential for their predictability. We employ a rolling window 
approach in order to capture the correlation dynamics for different window 
lengths and analyze the distributions of periods with statistically significant 
correlations. We find that for sufficiently large window lengths these distri- 
butions fit well to power-law behavior. We also measure the predictability 
itself by a hit rate, i.e. the rate of consistency between the signs of the ac- 
tual returns and their predictions, obtained from a simple correlation-based 
predictor. It is found that during these relatively brief periods the returns 
are predictable to a certain degree and the predictability depends on the 
selection of the window length. 

Keywords: Financial time series, episodic non-linearity, bicorrelation, time 
series prediction 

1. Introduction 

Recent empirical studies suggested that returns of financial time series 
depart from the random walk hypothesis by showing a certain degree of 
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long-term or short-term dependent relationships, thus violating the weak- 
form efficient market hypothesis [H, 0, 0] • The efficiency/non-efficiency can, 
for example, be assessed by evaluation of the Hurst exponent [3] or the ap- 
proximate entropy Q , as measures of long-term memory and randomness in 
time series. Since the efficient market hypothesis in its weakest form implies 
that the returns should be serially uncorrelated, its validity can also be as- 
sessed by looking for evidence of significant linear jo), 0] and non-linear Js|] 
serial autocorrelations (henceforth, correlations) in the returns. A systematic 
review of literature on the weak-form market efficiency of stock markets us- 
ing a wide array of statistical tests, such as the linear serial correlations, unit 
root, low-dimensional chaos, nonlinear serial dependence and long memory 
was provided by Lim and Brooks However, the market efficiency vary 
with time and the character of possible serial dependencies is not known. 
Moreover, they are of transient nature, showing up in random intervals just 
for a short time, which makes it difficult to exploit them for prediction pur- 
poses 

the non-linearities were inves 



In several previous studies [10|, [III, 1121, [13 



tigated by calculation of the bicorrelation test statistic, due to Hinich 
and performing the windowed-test procedure 
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15l . |16| . In order to capture 



the time variation, the data were split into a set of non-overlapping windows, 
the length of which was set to some ad-hoc fixed value. However, consider- 
ing the episodic and transient nature of the correlations, these may or may 
not be detected, depending on the window length used. Furthermore, the 
window length also influences the onset and offset of the significant correla- 
tions. Shorter lengths facilitate quicker response to changes in the correlation 
strength and can help pinpoint the arrival and disappearance of the tran- 
sient dependences, but on the other hand, they may lack adequate statistical 
power. 

In order to track the evolution of market efficiency over time, a rolling- 
window approach has been adopted to calculate a time-varying Hurst expo- 



nent 



171 . Il8l . Il9l |20| . Lim [21J has proposed a similar rolling window approach 



It 



for computing the bicorrelation test for various stock markets [22J, |23 
was shown that the market efficiency is not a static property that remains 
unchanged throughout the entire estimation period. Thus, examination of 
the presence of significant correlations in the rolling sample framework can 
provide a useful tool for tracking the changing degree of weak-form marker 
efficiency over time as well as for its ranking among different markets. In 
the currency exchange market, Brooks and Hinich 24J found that Sterling 
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exchange rates are characterized by transient epochs of dependencies sur- 
rounded by long periods of white noise. 

The objective of this paper is to study the dynamics of the linear and 
non-linear serial dependencies in financial time series, using a portmanteau 
test procedure in a rolling window framework. In particular, we focus on 
detection of episodes of statistically significant two- and three-point corre- 
lations in the returns of several leading currency exchange rates that could 
offer some potential for their predictability. We employ a rolling window 
approach in order to capture the correlation dynamics for different window 
lengths. In comparison with similar previous studies, our new contribution 
is using tools from statistical physics in order to analyze distributions of 
the periods with statistically significant correlations and relate them to the 
predictability obtained from a simple correlation-based predictor. 



2. Data 

We use hourly average exchange rates per Euro (EUR) on the US dollar 
(USD), Canadian dollar (CAD), Swiss franc (CHF), British pound (GBP) 
and Japanese yen (JPY). The sample period is from March 1, 2004 to De- 
cember 17, 2009 (a total of 37,620 observations). In Fig. [TJ we plot the 
nominal exchange rates of the respective currencies vs the Euro. In Figs. [2] 
and [3] we show first logarithmic differences of the nominal exchange rates and 
their frequency distributions, respectively. The histograms and the summary 
statistics, presented in Table [TJ indicate that the data deviate from the nor- 



mal distribution. This observation is also supported by the Jarque-Bera [25 
normality tests, which reject the null hypothesis of normality with very low 
p- values. 



3. Methods 

3.1. Detection of significant correlations and bicorrelations in a rolling win- 
dow 

Let R(t) represent the logarithmic returns of the length N, defined by 
R(t) = \n(P(t + 1)/P(t)), where P(t) are the nominal exchange rates and t 
is a time label, t = 1,...,N. Within the rolling window framework we test 
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-0.0261 


Std deviation 


1.3e-03 


1.4e-03 


0.7e-03 


l.le-03 


1.7e-03 


Skewness 


-0.1521 


-0.0486 


-0.3913 


-0.3651 


0.2721 


Kurtosis 


17.59 


13.14 


52.95 


22.03 


27.41 


JB normality test statistic 


333,840 


161,300 


3,912,000 


568,380 


934,710 


p- value 


< 1.0e-03 


< 1.0e-03 


< 1.0e-03 


< 1.0e-03 


< 1.0e-03 



Table 1: Summary statistics of logarithmic returns 

for both linear and non-linear correlations in a window of a specified length 
n, that is rolling one point forward eliminating the first observation and 
including the next one, until the last observation. The data in each window 
are standardized to have zero mean and unit standard deviation, as follows: 

OR 

where Mr and Sr are the sample mean and standard deviation of the win- 



dow. The window-test procedure |15|, [16J uses correlation and bicorrelation 



potmanteu tests for the detection of linear and non-linear serial dependen- 
cies within a time window. The null hypothesis is that the standardized data 
x(t) in each window are realizations of a stationary pure white noise process. 
The alternative hypothesis is that the process has some non-zero correlations 
C xx {r) = E[x(t)x(t + r)] or bicorrelations C xxx (r, s) = E[x(t)x(t + r)x(t + s)] 
within 0<r<s<L<n, where L is the number of lags. The sample 
correlations and bicorrelations in a rolling window are calculated as follows: 

n—r 

C xx (r) = - -5>(t)x(t + r), (2) 

j n—s 

C xxx (r, s) = } x(t)x{t + r)x(t + s). (3) 

v > t=i 
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The test statistics in a rolling window for non-zero correlations H xx and 
bicorrelations H xxx , are respectively given by: 

L 

lf« = £(n - r)<£(r) ~ ^(L) (4) 



and 



L s-1 



H xxx = s)C 2 xxx (r, s) ~ X 2 ((L - 1) V2). (5) 

s=2 r=l 

We consider several window lengths n = 2 k , k = 3, ...,10. The number of 
lags is specified as L = n b , where < b < 0.5 is a user-specified parameter. 
In order to account only for the presumably most relevant lags, we chose a 
fixed value of L = 2 (i.e. r — 1, s — 2), which, for example, for n = 2 4 



corresponds to b = 0.25, as used in Ref. [10 . 

We note that in order to identify non-linearities by the H xxx test statistics, 
data pre-whitening is necessary. This is achieved by filtration of the linear 
component by an autoregressive AR(p) fit. The p order of the AR(p) model 
is chosen between 1 and n — 1 as the optimizer of Schwarz's Bayesian Crite- 



rion 



26|, |27j. The null hypothesis of linear /non-linear correlation is accepted 



or rejected in each window at a risk level of 5%. 

3.2. Simple correlation-based predictor 

Above we identify episodes with significant correlations which indicate 
transient inefficiency of the market. Emergence of such episodes offers some 
opportunity for the prediction of the price change. In the following we define a 
simple predictor of the price change direction in the intervals with significant 
correlations. Prediction of the price change direction I (+1 - increase, —1 - 
decrease, - no prediction) at the time step t + 1 is based on the assumption 
of consistency between the signs of the neighboring price changes and the 
sample correlation in a window: 

I(t + 1) = sigR(C xx (l)x(t))6(a - Pxx ), (6) 

where p xx is a p-value of the H xx statistic, 6 is the step function {6 = 1 
if p xx < a and otherwise), and a determines the confidence level for no 
correlation null hypothesis rejection (we used a = 0.05). Hence, if the sample 
correlation in a given time window with significant correlation is positive 
(negative) the predictor predicts the same (opposite) price change direction 
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as the last observed change. The predictor's performance is evaluated by a 
hit rate, i.e. the rate of consistency between the signs of the actual returns 
and their predictions. We note that we avoided using a more complicated 
prediction model that would involve parameters. Instead, we tried to link the 
presence of significant correlations with the degree of predictability as directly 
as possible, thus eliminating the influence of the parameters' inference on 
the prediction. By using the current simple model the only relevant issue 
is weather the magnitude of the estimated correlation (positive or negative) 
significantly differs from zero or not. 



4. Results 

4-1. Frequency and cluster size distribution of significant H xx and H xxx win- 
dows 

The arrival of periods with significant correlations and their lengths are 
stochastic processes and their evaluation can shed some light on the overall 
efficiency of a particular market as well as the interpretation of the events, 
such as crises, government regulations, etc., that could have lead to the in- 



crease or decrease of the efficiency [2Cj, [22J, |23| . Lim et al. |22|, |23| proposed 
the percentage of time windows in which the market departs from efficiency 
as an indicator for assessing the relative efficiency. In Fig. H]we plot the per- 
centage of the significant (a) H xx and (b) H xxx rolling windows for different 
currency pairs as a function of the window length n. We can observe that 
the percentage of the significant H xx and H xxx windows tends to increase 
with increasing n, with some currency pairs (e.g., CHF/EUR) being more 
and some less (e.g., USD/EUR) sensitive to the variation of n. Based on the 
efficiency indicator proposed by Lim et al. [221, 23[, the largest (smallest) ef- 



ficiency is found in USD/EUR (CHF/EUR) pair, featuring up to 13.5 (22.2) 
% of the significant H xxx windows for the considered values of n. 

The epochs of the relative inefficiency represent clusters of instants with 



significant correlations of various sizes. As pointed out in [10l . |24j , the de- 
pendency structures are short-lived and not persistent enough to be properly 
exploited for prediction purposes. Therefore, in this perspective it is in- 
teresting to look into the size distribution of the clusters with significant 
correlations. In Fig. [5] cumulative distributions are plotted in a log-log scale 
for the CAD/EUR pair with different values of the window length n. We 
can see qualitatively similar behavior in both the significant H xx and H xxx 
cluster size distributions. Namely, for sufficiently large n, fat tails can be 
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p- value 




p- value 


64 


4 1.77 


0.00 


19 2.82 


0.01 


128 


13 1.97 


0.03 


19 2.31 


0.08 


256 


12 1.8 


0.12 


17 1.95 


0.14 



Table 2: Power-law fits to distributions of clusters of significant H xx and H xxx windows 
for the CAD/EUR pair for different window lengths n, using the maximum- likelihood 
fitting method with goodness-of-fit tests based on the Kolmogorov-Smirnov statistic [28| . 
The parameters x m i n and a correspond to the estimates of the lower-bound of the power- 
law behavior and the scaling exponent, respectively. Statistically significant values with 
p > 0.1 arc in bold. 



observed. As evidenced from Table [21 a reasonable fito to the power-law be 
havior, with the estimated exponent a ~ 1.8 for H xx and a ~ 1.95 for H. 



XXX ■ 



can be obtained by the maximum likelihood method 28] for n = 256. This 
behavior is in clear contrast with the cluster-size distributions of Gaussian 
white noise, depicted in Fig. [5] by the square symbols. Similar distributions 
can also be observed using the remaining currency pairs, as shown in Fig. |6j 

4-2. Predictability in significant H xx windows 

The epochs with significant correlations can be used to forecast the next 
price change by using the correlation-based predictor ([6]). Its performance 
in relation with the significance of the linear correlations is demonstrated 
in Fig. [7J for a selected period of CAD/EUR data. The cumulative hit 
rate displays a sharp increase during the epochs characterized by persistent 
contiguous clusters with very low p- values of the H xx statistic (e.g., the epoch 
within the green bars). The hit rate typically deteriorates at the end of 
these epochs due to the fact that, especially for larger window lengths n, 
the response to changes in the correlation strength is relatively slow and 
the low p- values persist for some time even after the transient dependences 
have disappeared. A quicker response can be achieved using smaller window 
length n, but then the epochs with transient dependences tend to become 



1 Normally p- value indicates that the null hypothesis is unlikely to be correct and there- 
fore low values are considered good. In [2^], by contrast, the p- value is used as a measure 
of the hypothesis that is being verified, and hence high values are good. The power law is 
a plausible hypothesis for the data if the resulting p- value is greater than 0.1. The number 
of semi-parametric bootstrap repetitions of the fitting procedure is set to 1000. 
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shorter, as already evidenced in Fig. [5J On the other hand, the predictions 
made in the epochs with p-values below the threshold a = 0.05 but with 
higher variability are not as successful (e.g., the epoch within the red bars). 
The tendency of obtaining a better hit rate as the cluster size increases is 
demonstrated in Fig. 8(a) for the CAD/EUR pair and different values of the 
window size n. An excellent hit rate is achieved in relatively large clusters 
for smaller values, such as n = 8, 16. However, as shown in Fig. El such 
clusters are scarce and the overall hit rate is governed by the dominant small 
clusters that show poorer hit rate. For larger values of n (e.g., n = 256), the 
fat tails tell us that large clusters are more ubiquitous, nevertheless, the hit 
rate in these large clusters is markedly worse than in the those produced by 
small n. As a result, looking at the total hit rate, no significant differences 
are observed between different values of the window size n. In Fig. 8(b) , we 



again show qualitative difference compared to the uncorrelated random time 
series, which does not not feature larger clusters with significant correlations 
and the hit rate is not affected by the cluster size but randomly fluctuates 
around the value of 0.5. 



5. Conclusions 

The objective of this paper was to study the dynamics of the linear and 
non-linear serial dependencies in financial time series, using a windowed port- 
manteau test procedure in a rolling window framework. In particular, we fo- 
cused on detection of episodes of statistically significant two- and three-point 
correlations in returns of several leading currency exchange rates that could 
offer some potential for their predictability. Since the correlations, which 
can be viewed as an indicator of market efficiency, are not static but they 
evolve in time, we employed a rolling window approach in order to capture 
their dynamics for different window lengths. By analyzing the distributions 
of the periods with statistically significant correlations it was found that the 
percentage of the periods with statistically significant two- and three-point 
correlations in a window of past n observations varies with n and differ- 
ent currencies show different degrees of the dependences. For example, the 
CHF/EUR pair shows a relatively high percentage of significant correlations, 
which increases with the window length, while the USD/EUR pair displays 
a relatively small percentage with little sensitivity to the window length. We 
found that for sufficiently large window lengths these distributions of the pe- 
riods with statistically significant correlations fit well to power-law behavior. 
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We also measured the predictability itself by a hit rate, i.e. the rate of consis- 
tency between the signs of the actual returns and their predictions, obtained 
from a simple correlation-based predictor. It was found that during these rel- 
atively brief periods the returns are predictable to a certain degree and the 
presence of large contiguous clusters of instants with significant correlations 
increase potential for predictability. 

In the present study we did not attempt to use the detected non-linearities 
in the windows with significant bicorrelations for prediction purposes. Pre- 
diction of such non-linear processes is a challenging task, nevertheless, some 
attempt using a simple model has already been made [29|. In a view of 
the above, it would be interesting to employ jointly linear and non-linear 
prediction models during the periods featuring solely linear (non-zero corre- 
lation and zero bicorrelation) and non-linear (zero correlation and non-zero 
bicorrelation) dependencies, respectively. 

References 

[1] A. W. Lo, Econometrica 59, 1279 (1991) 

[2] C. Eom, S. Choi, G. Oh, W. S. Yung, Physica A 387, 4630 (2008) 

[3] C. Eom, G. Oh, W. S. Yung, Physica A 387, 5511 (2008) 

[4] H. E. Hurst, Trans. Am. Soc. Civil Eng. 116, 770 (1951) 

[5] S. M. Pincus, Proc. Natl. Acad. Sci. USA 88, 2297 (1991) 

[6] E. F. Fama, Journal of Business 38, 34 (1965) 

[7] C. W. J. Granger, O. Morgenstern, Kyklos 16, 1 (1963) 

[8] M. J. Hinich, D. M. Patterson, Journal of Business and Economic Statis- 
tics 3, 69 (1985) 

[9] K. P. Lim, R. D. Brooks, Journal of Economic Surveys 25, 69 (2011) 

[10] A. Serletis, A. G. Malliaris, M. J. Hinich, P. Gogas, Open Econ Rev, 
DOI:10.1007/sll079-010-9194-9 

[11] C. Brooks, M. J. Hinich, J. Empir Finance 20, 385 (1999) 



9 



[12] C. A. Bonilla, R. Romero-Meza, M. J. Hinich, Applied Economics Let- 
ters 13, 195 (2006) 

[13] C. A. Bonilla, R. Romero-Meza, M. J. Hinich, Applied Economics 39, 
2529 (2007) 

[14] M. J. Hinich, J Nonparametric Statist 6, 205 (1996) 

[15] M. J. Hinich, D. M. Patterson, Mimeo, University of Texas at Austin 
(1995) 

[16] M. Hinich, D. M. Patterson, In: M. T. Belongia and J. M. Binner (Eds.), 
Money, Measurement and Computation, (Palgrave Macmillan, London, 
2005) 61 

[17] D. O. Cajueiro, B. M. Tabak, Chaos, Solitons & Fractals 22, 349 (2004) 

[18] D. O. Cajueiro, B. M. Tabak, Physica A 336, 521 (2004) 

[19] J. Alvarez-Ramirez, J. Alvarez, E. Rodriguez, G. Fernandez-Anaya, 
Physica A 387, 6159 (2008) 

[20] Y. Wang, L. Liu, R. Gu, J. Cao, H. Wang, Physica A 389, 1635 (2010) 

[21] K. P. Lim, Physica A 376, 445 (2007) 

[22] K. P. Lim, R. D. Brooks, J. H. Kim, International Review of Financial 
Analysis 17, 571 (2008) 

[23] K. P. Lim, R. D. Brooks, Chaos, Solitons and Fractals 40, 1271 (2009) 

[24] C. Brooks, M. J. Hinich, Applied Economics Letters 5, 719 (1998) 

[25] C. M. Jarque, A. K. Bera, Econ Lett 6, 255 (1980) 

[26] H. Lutkepohl, Introduction to Multiple Time Series Analysis, 2nd edi- 
tion, (Springer- Verlag, Berlin, 1993) 

[27] A. Neumaier, T. Schneider, ACM Transactions on Mathematical Soft- 
ware 27, 58 (2001) 

[28] A. Clauset, C. R. Shalizi, M. E. J. Newman, SIAM Review 51, 661 
(2009) 

10 



[29] D. Sornette, V. F. Pisarenko, Physica D: Nonlinear Phenomena 237, 429 
(2008) 



11 




x 10 



(a) USD/EUR 



0.7 
0.68 
0.66 
0.64 
0.62 

0.6 
0.58 



\ 



1 



x 10 



(c) CHF/EUR 

x 10" ; 



9 
8.5 

8 
7.5 

7 
6.5 

6 

5.5, 



HfK 



1 



0.75 



0.7 



0.65 



0.6 



0.55. 



1.6 



x 10 



(b) CAD/EUR 



1.4 
1.3 
1.2 
1.1 

1 




x 10 



(d) GBP/EUR 



x 10 



(e) JPY/EUR 



Figure 1: Nominal exchange rates. 



12 



0.02 




x 10 



(a) USD/EUR 



0.02 



0.01 



-0.01 



-0.02. 




x 10 



(b) CAD/EUR 



0.02 



0.01 



-0.01 



-0.02 




x 10 



0.02 



0.01 



-0.01 



-0.02„ L 




x 10 



(c) CHF/EUR 



(d) GBP/EUR 




x 10 



(e) JPY/EUR 



Figure 2: Logarithmic returns. 
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Figure 3: Frequency distributions of the logarithmic returns with the Gaussian curve 
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Figure 4: Percentage of the significant (a) H xx and (b) H xxx rolling windows for different 
currency pairs as a function of the window length n. 
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Figure 5: Distributions of clusters of significant (a) H xx and (b) H xxx windows of size s, 
for the CAD/EUR pair for different values of the window length n. The power- law (PL) 
fit is done for n = 256, using the maximum-likelihood method (2||. The filled squares 
correspond to a simulated random time series ~ N(0, 1) and different window lengths. 
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Figure 6: Distributions of significant (a) H xx and (b) H xxx cluster sizes s for different 
currency pairs for the window length n — 256. 
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Figure 7: Evolution of (a) the p- value of the H xx statistic (p xx ) and (b) the hit rate in 
the periods of significant correlations. The data correspond to the first 18,000 CAN/EUR 
return samples and the window size n = 256. 
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Figure 8: Prediction hit rate as a function of the cluster size s for (a) the CAD/EUR pair 
and (b) a simulated random time series ~ N(0, 1), using different window lengths n. 
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